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■ In reactions with qq production or qq annihilation, initial- and final- state 

' interactions give rise to large corrections to the lowest-order cross sections. We 

Q_i' evaluate the correction factor first for low relative kinetic energies by studying 

Oh: 



the distortion of the relative wave function. We then follow the procedure of 
Schwinger to interpolate this result with the well-known perturbative QCD 



^ ■ vertex correction factors at high energies, to obtain an explicit semi-empirical 



correction factor applicable to the whole range of energies. The correction 
factor predicts an enhancement for qq in color-singlet states and a suppression 
for color-octet states, the effect increasing as the relative velocity decreases. 
Consequences on dilepton production in the quark-gluon plasma, the Drell- 
Yan process, and heavy quark production processes are discussed. 
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I. INTRODUCTION 



The possibility of nuclear matter going through a phase transition into a deconfined 
quark gluon plasma (QGP) state during high-energy heavy-ion collisions makes these colli- 
sions the focus of intense experimental and theoretical research |l|]. In these collisions, the 
constituents of the possible QGP or the partons of the colliding nucleons can react at varying 
energies. Their reaction products such as dileptons and photons provide the signals and the 
backgrounds for the detection of the quark-gluon plasma. The magnitudes of the signals 
depend on the cross sections for dilepton and photon production from the constituents of 
the plasma. The rates of approach to chemical equilibrium and thermal equilibrium of the 
plasma depend also on the cross sections for reactions among the constituents. 

When quarks, antiquarks and gluons interact, they interact as constituents in the quark- 
gluon plasma or as partons in the colliding nucleon. However, the lowest-order Feynman 
diagrams and the next-to-leading order Feynman diagrams are the same for the basic reaction 
processes involving g, q and gluons, whether the reaction takes place in the environment of 
the quark-gluon plasma or in the environment of partons in nucleon-nucleon collisions. In 
these basic reaction of quarks, antiquarks and gluons, the next-to leading order diagrams, 
including the initial- and final-state interactions and gluon radiations, give rise to large 
corrections to the lowest-order cross sections. It will be useful to develop an analytical semi- 
empirical correction factor for the basic reaction processes over the whole range of relative 
energies so that in the next level of approximation the basic lowest-order reaction cross 
sections of quarks, antiquarks and gluons can be corrected on the same footing. Additional 
effects and refinements such as the plasma screening and temperature can be added on in 
the future as our theoretical understanding is developed further. 

For processes involving qq production and annihilation at high energies, it has been al- 
ready generally recognized in perturbative QCD that the lowest-order Feynman diagrams 
give only an approximate description ||23|. In the case of the Drell-Yan process for 
instance, the experimental cross section attributed to the process is about a factor of 2 
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to 3 greater than what one predicts by using the lowest-order Feynman diagram. A phe- 
nomenological i^'-factor is introduced by which the lowest-order results must be multiplied in 
order to bring the lowest-order QCD predictions into agreement with experiment. Higher- 



order QCD corrections to the Drell-Yan cross section, up to 0{as) [|TT| and 0{al) [jT3 
have been worked out. These investigations show that in the high energy limit of massless 
quarks, the K factor can be well accounted for by including QCD corrections [p!Ii [l3|. The 
most important contribution to the dilepton cross section due to the Os-order (the next-to- 
leading-order) diagrams is the vertex correction at the qq'y* vertex involving the initial-state 
interaction between q and q and the gluon radiation from q and q. According to ||8|-[10|, the 



vertex correction leads to a correction factor equal to (1 + 27ras/3). For a strong interaction 
coupling constant = 0.3, this vertex correction gives a factor of about 1.7. Additional 
contributions from the Og-order Compton diagrams bring the Drell-Yan i^-factor within the 
observed range of 1.6 to 2.8. For charm and heavy-quark production, the cross section for 



the lowest-order QCD was given in [|T^,|T^. Higher-order QCD corrections to the cross sec- 
tions and single-particle inclusive differential cross sections have also been obtained P^7|-pD|. 
If one includes the next-to-leading-order of QCD, the i^-factor is found to range from 2 to 
3, depending on which heavy flavor is produced. The threshold behavior for heavy-quark 
production due to the distortion of the Q-Q relative wave function has been examined in 
Refs. |@-|^. The correction is quite large near the threshold. 

We focus in this work on basic reaction processes involving q and q in the initial and 
final states. Other reactions involving gluons are important part of the dynamics in the 
plasma or partons and will be the subject of our subsequent studies. We investigate here 
the corrections to the tree-level cross sections involving the production and annihilation of 
qq at all relative energies. For low energies, we seek an analytical vertex correction factor 
arising from wave function distortion by virtue of the q-q color potential. The distortion 
correction depends on as/v where v is the magnitude of the asymptotic relative velocity. It is 
the most important correction to the tree-level descriptions at energies where masses cannot 
be ignored, as has been earlier investigated for heavy quark production processes At 



high energies, we use the well-known PQCD results which include the initial- or final-state 
interactions in addition to real gluon radiations. We shall seek an interpolation to join the 
vertex correction factor from low energies to these well-known perturbative QCD vertex 
correction factors at high energies. For the purpose of interpolation, we shall include the 
proper relativistic kinematics and follow the interpolation procedure suggested by Schwinger 
|30| . We obtain a simple analytical semi-empirical correction factor which can be used over 
the whole range of energies for both the annihilation and the production of a quark-antiquark 
pair. 

This paper is organized as follows. In Section II, we first give a general discussion on the 
correction factor at low relative velocities arising from the distortion of the wave function 
at the point of qq annihilation or production. This correction factor is found to be an 
enhancement factor for the color-singlet states and a suppression factor for the color-octet 
states. The interpolation of this result with the known perturbative QCD correction factor 
at high energies provides correction factors which can be applied over the whole range of 
relative velocities. In Section 3, we show how we can use vertex correction factors to improve 
the lowest-order Drell-Yan cross sections. We discuss the use of these correction factors for 
dilepton production in the quark-gluon plasma in Section 4. The correction factors lead 
to an enhancement of the dilepton production probability in the plasma. In Section 5, we 
examine the qq —>■ QQ process, where both qq annihilation and QQ production are found 
to be associated with suppressive correction factors. Section 6 summarizes the present 
discussions. 

II. GENERAL CONSIDERATION OF THE CORRECTIVE FACTOR FOR qq 

ANNIHILATION AND PRODUCTION 

What is the effective potential between the quark and the antiquark for continuum states 
that populate the production and annihilation vertices? The processes of production and 
annihilation are characterized by a region of very small relative q-q separations, at a linear 
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distance of ~ a/y/s for a virtual photon intermediate state and a distance of ~ as/y/s for 
a virtual gluon intermediate state. The effective part of the q-q interaction in this region 
of small relative separation associated with the annihilation or production process is the 
inverse-r, Coulomb-like term. The linear part of the potential, that serves to effect the 
binding and impose confinement, controls the large r behavior and needs not be considered 
in the collapsed spatial zone comprising the annihilation or production vertex. Therefore, we 
consider the interaction between a quark and an antiquark with an invariant mass Q = y/s 
as described by an effective Coulomb-type potential A = {Aq, 0) involving their relative 
coordinate r f^-^ 

Aoir) = (2.1) 
r 

where agfr is the effective strong- interaction coupling constant, related to the strong inter- 
action coupling constant as by the color factor Cf, 

Oeff = Cftts . (2.2) 

The running coupling constant is 

127r 

~ (33 - 2nf) ln(gVA2) ' ^ ' 

where the flavor number Uf can be parametrized to be 

q=u,d,s,c,b,t 



with the step function 6 and the threshold behavior inferred from that of qq QQ cross 
section. 

When q and q are in their relative color-singlet states, the interaction is attractive and 
the color factor Cf is 4/3 When q and q are in their relative color-octet states, the 

interaction is repulsive and the color factor Cf is —1/6 [Q. 

It is worth mentioning that the notion of a potential is a useful concept, as one can 
infer from the success of the quasi-potential of Tordorov l2^ and Crater and van Alstine's 



relativistic constraint dynamics to describe qq and e^e as a relativistic two-body problem 
26| , |27[] . In these descriptions of constraint dynamics, the compatibihty of the two-body 



equations of motion requires that the two-body potential must be expressed in terms of 
the relative four-coordinate orthogonal to the total four- momentum P. It is therefore 
convenient to work in the center-of-mass coordinate system, where P = {^/s, 0) and the 
four-coordinate x± is the relative coordinate r used in Eq. (|2.1| ). In this coordinate system, 



the two-body system is described by an energy for relative motion given by 

„ s - 2ml 



and the relativistic generalization of the reduced mass of the system is = / y/s. The 
corresponding value of the on-shell relative momentum squared at r ^ oo is 

p' = b\s) = el-ml='^^^. (2.5) 



(See Eqs. (2.13a), (2.13b), and (2.13c) of Ref. 

A proper treatment of the relativistic two-body problem will require the solution of the 
two-body Dirac equation as formulated in Ref. [ PB| , P7[ which is rather complicated. To gain a 
simple insight into the effect of the interaction between the quark and the antiquark, we shall 
treat the problem approximately as a single fermion spinor in the external field of Eq. (p.l|) 



for which the solution is already known p4[, with the stipulation of the correct relativistic 
kinematics to describe the relative motion between the quark and the antiquark. As is well- 
known, upon using the usual approximation as given by Ref. |2^, the wave function which 



satisfies the correct boundary condition of an incident plane wave is 

^ = A^e*P-''(l - —a ■ V)uF{i^, 1, i{\p\r -p-r)), (2.6) 

where u is the spinor for a quark, F is the confluent hypergeometrical function, N is the 
normalization constant given by 

1^1' = (2-7) 
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and the parameter ^ (which can be positive or negative) is related to the couphng constant 
tteff and the magnitude of the asymptotic relative velocity v by 

e = — • (2.8) 
The magnitude of the asymptotic relative velocity v is the ratio of the momentum b to the 



energy in the system of relative coordinates. Following Todorov |25| and Eqs. (21.13a) 



2.13c) of Crater et al ||2^, the relative velocity for the quark and the antiquark in their 



center-of-mass system obtained by using Eqs. (|2.4|) and (|2.5|) is 



^ '^^ (2.9) 



s — 2m^ 



which gives v ~ 2^1 — Am'^/s when ~ 2mq and v 1 when s — > oo. The quantity 
V is the asymptotic relative velocity obtained with the proper relativistic kinematics and 



differs from the quantity (3 = y 1 — 4m^/s defined by Fadin et al. 0. For the Drell-Yan 
qq annihilation processes or for the heavy-quark pair production processes, a quark and an 
antiquark are first annihilated into (or produced from) a virtual photon or a virtual gluon. 
As annihilation (or production) takes place in the vicinity of the region around r = 0, the 
probability for the annihilation process is proportional to the absolute square of the wave 
function at r = 0. When one averages over the spins of the quark, the square of the wave 
function at the origin in the case with interaction Aq is 



\m\io = N' 



, /acfr\^ 2(1 - cos6') 

' + l:^J^ — 2 — 



(2.10) 



where 9 is the angle between the incident momentum p and r. After averaging over 9, we 
have 

WO)li = ^^^(l + %. (2.11) 

The above result was obtained by treating the qq system approximately as a fermion spinor 
in the field ( |2.1D of a spinless particle and the correction term tteff/2 inside the bracket arises 
from the spin of one of the two particles. When the spin of the other particle is taken into 



account, there is an additional «ofj/2 contribution, and the square of the wave function at 
the origin is modified to 

1^(0)1^0 = ^^^(l + «e\). (2.12) 

On the other hand, in the perturbative expansion in which there is no interaction between the 
quark and the antiquark in the lowest-order approximation, the square of the wave function 
|'0(O)|q is unity. Therefore, by considering the interaction between the quark and the anti- 
quark as arising from an effective Coulomb-like interaction, the annihilation or production 
cross section is modified by the factor, which is a generalization of the 'Gamow-Sommerfeld' 
factor [^,^, given by 



|^(0)|i„ 27re , , 

The expansion of the above K factor arising from the distorted wave function gives 



K 



^ 2< (27rO^ (27rO^ (20^ (20^ 
2 12 720 30240 1209600 



l + O, (2.14) 



which contains higher order terms with alternating signs. It is clear that for the case when 
\TiCfas/v\ is much greater than 1, the convergence of the perturbation expansion as a power 
of Us will be very slow. Many higher-order terms are needed to approach the full result Eq. 
( p.l3| ). For that case, a perturbative treatment of the correction factor is not useful and 
a non-perturbative treatment such as presented here in terms of wave-function distortion 
represented by Eq. ( |2.13|) is a meaningful concept. 

The correction factor Eq. (|2.13|) has been obtained by considering only the exchange of 
virtual gluons between the quark and the antiquark. The emission of real soft gluons is not 
taken into account. Because radiative corrections involving the emission of real gluons is 
unimportant for qq systems with small relative velocities, the results of Eq. ( |2.13| ), without 
the inclusion of real soft gluons, is applicable to the region of low relative velocities. 

At very-high relative velocities, besides the distortion effect arising from the initial-state 
and final-state interactions between the quark and the antiquark, it is necessary to include 
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the radiation of real gluons. In that case, the amphtude from the emission of soft gluons 
interferes with the amphtude arising from the exchange of a virtual gluon between the quark 
and the antiquark. At very high relative velocities for which the masses of the quarks can be 
neglected, perturbative QCD calculations have been carried out by including the exchange of 
a virtual gluon between the quark and the antiquark and the emission of real gluons from the 
two particles [|r2,|l^,|,^ . These high-energy vertex correction i^-factors are different for qq 
annihilation or for qq production. As the vertex correction gives the dominant contribution, 
we shall consider the approximation where the phenomenological correction factor K is given 
just by the vertex correction. Up to the first order in as in the high-energy limit of massless 



quarks, the vertex i^-factor for qq annihilation is JlT 

K= 1 + vraeff \7r^ + a 
= l + vraeff X 0.884, (2.15) 

which is about 2 for a typical Oefr = 0.4 for color-singlet states. The vertex correction 
i^'-factor for qq production in the e+e" qq process is [^,^ 

K = l + — , (2.16) 

71 

which is about 1.1 for a^s = 0.4. On account of the fact that the singlet color factor for 
e~^e~ annihilation is 4/3, the correction factor for qq production is 

K = l + ^. (2.17) 

Thus, for qq in a color-singlet state at high energies, the high-order QCD corrections are 
large for qq annihilation and small for qq production. In contrast, at low relative velocities, 
Eq. ( |2.13| ) gives equally large QCD corrections for qq annihilation or production. 



Knowing the two limits of the correction factors from Eqs. ( |2.13| ), ( |2.15| ), and ( p. 171 ) 



we can follow the interpolation procedure suggested by Schwinger [^] and used by Barnett 
et al. 0]. We can obtain a semi-empirical correction factor which can be used for the whole 
range of relative velocities. For qq annihilation in the initial state, we define a function 
/W(^) by 



eff 



1/15 
- + v(-l + — + - 



(2.18) 



For qq production in the final state, we define a function f^^\v) by 



eff 



- + v{-l + — 



(2.19) 



Associated with the annihilation or production of a qq pair at a center-of-mass energy ^/s, 
we propose a semi-empirical vertex K-factor as 

27r/(*'-^)(t;) 



K(''f\q) 



(2.20) 



1 - exp{-27r/(*'/)(t;)} 

where the fiavor label q in K^^\q) is included to indicate that K^^'^ depends on the quark 
mass m^, the superscript {€) and (/) denotes qq initial-state annihilation and /inal-state 
production respectively, and the magnitude of the relative velocity v is given by Eq. ( p.9| ). 
This vertex correction factor agrees with the results of Eq. ( |2.13| ) for low relative velocities. 
For high relative velocities and |vrae//| << 1, Eq. ( |2.20D gives 



K^iJ) 1 + tt/ 



(2.21) 



and Eq. ( ^:20D agrees with the well-known results of Eqs. ( |2.15|) and ( |2.17D from perturbative 
QCD for high relative velocities, up to the first order in Thus, Eq. ( p.20|) can provide 
a good description of the vertex corrections for qq systems over the whole range of relative 
velocities, for qq annihilation or qq production. Because the vertex correction gives the dom- 
inant QCD correction, the vertex correction factor ( |2.20| ) can be used to give the corrections 
to lowest-order QCD results when a qq pair is annihilated or produced. It could equally 
well be used with a subscript 'v' to identify its association with the vertex correction and 
distinguish it from the usual phenomenological K factor that usually quantifies the QCD 
corrections in the massless limit only. However, it has been retained as K for compactness 
and to allow easy applicability of the correction. 
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III. DRELL-YAN PROCESSES 



As an example of the application of the correction factor K^'^'-^^ of Eq. ( p.2(]| ), we can use 
it to improve theoretical tree-level estimates of the Drell-Yan cross section. In the Drell-Yan 
process, a qq pair is annihilated in the initial state, and the appropriate correction factor 
is K^'^\ As the virtual photon in the qq'j* vertex selects the color-singlet combination of 
the annihilating qq pair, the corresponding color factor C/ is 4/3, and the effective coupling 
constant needed to calculate K^^^ in Eq. (13) is aefr = Aos/S. 

If only the lowest theoretical estimate da^'^ / dQ'^dy is available, an improved determina- 
tion of the Drell-Yan cross section for qq — > can be provided by 



dQ'^dy dQ'^dy 

In Fig. 1, we show K^^\q) as a function of the invariant mass Q = y/s of the dilepton pair 
calculated with A = 0.3 GeV. For a qq pair in a color-singlet state at high relative velocities, 
the K^^^ factor is about 2.2 which agrees with the experimental value of 1.6 to 2.8 in the 
Drell-Yan process. For color-singlet states with low relative velocities, the correction factor 
is quite large. The correction factor for a qq pair in the color-octet state approaches 0.9 as 
y/s increases and it deviates from this constant behavior at low relative velocities when the 
energy is near the qq mass threshold. 

When the next-to-leading order results, including other diagrams such as Compton dia- 
grams are available, the estimate of the Drell-Yan cross section can be improved to include 
the distortion corrections as well. By following a procedure which was suggested by Harris 
and Brown ^l|], we can avoid double counting. 

Accordingly, the next-to-leading order results of Kubar et al. for the Drell-Yan cross 
section can be modified to include distortion effects by using the i^-factor as given by 
Eq. (|2.13|) prior to interpolation with the high energy QCD corrections. Using most of the 
notation of Kubar et al., the cross section for the production of an pair with an invariant 
mass squared and a rapidity y for the collision of two hadrons is 
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da da^^ da^ da^ 



with 



dQ'^dy dQ'^dy dQ'^dy dQ'^dy ' 



dcr^ r d(7^ 



(3.1) 



(3.2) 



where i = LO, A, C represent the contributions from the lowest order Drell-Yan diagram, 
the annihilation diagrams, and the Compton diagrams, respectively. The function is the 
product of the quark and antiquark structure functions of the colliding hadrons. 



Qf{tiM) = gfo (^1)^20(^2) + gio(^i)?2o(^2 



(3.3) 



and / is the flavor label. 

When the distortion due to the wave function is taken into account, the result of Kubar 



et al 0] can be modified to be 

da^o da^ 
dQ^dy dQ'^dy 



9Q2sV e 



6{h - Xi)S{t2 - X2)K^'^ 



X 1 + 



277 



X1X2 



laeffl 



6{t2 - X2] 



+ In 



+ 21n^^ln^^- 

(1 - a;i)(l - 0:2) 1-xi I-X2 

2xi(l - X2) 



\ — In 7 TT 

^ 2 l-a;i 1- 



+ 



tl{ti-Xi)+ X2{ti-Xi) 2{ti-Xi)+ ti tf 



2 3a;i\ 

---^) + (l-2) 



^^(^1,^2) 



+H^ih,t2: 



(3.4) 



TT \[iti- Xi){t2- X2)] , 

where K^^^ is now given by Eq. ( p.l3| ), the distributions l/(t — x) + , and l/[(ti — Xi)(t2 — 2^2)] + 
are defined by 



X [t — X). 



t — X 



and 



f{h,t2 



1 /•! 

XI ^ Jx2 ^ [(tl - Xi) {t2 - X2)] + 



^ ^+ /(^l' ^2) - /(^l, ^2) - f{XlM) + f{xi,X2) 
1 / '^2 77 TTT ^ 

xi Jx2 [J'l ~ Xi)\l2 — X2) 



In another perturbative treatment of the Drell-Yan process, an estimate of the first-order 
vertex correction leads to a vertex K factor of the form |^^10[ 
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K' =1 + 27ras/3 , (3.5) 

and the generalization to higher-orders of the form 

K' = e^™^/^ . (3.6) 

However, this result has been obtained only in the high energy limit and may not be applied 
to the cases of low dilepton invariant masses. Furthermore, the vertex correction factor of 
Eq. ( p.5|) appears to be only a part of the whole vertex correction factor of Eq. (|2.15|) . 

Combined with the running coupling constant, the use of the correction factor Eq. (|2.20|) 
allows us to probe lower dilepton invariant masses and the important low- a; region with 
better justification than massless limits. 



IV. THERMAL DILEPTON PRODUCTION 

The importance of the possibility of hadron matter going through a phase transition 
into a deconfined quark-gluon plasma state, during high-energy heavy-ion collisions has 
been highlighted earlier While the basic process leading to the formation of dileptons 
is the same as that for the Drell-Yan process, the distributions and the characteristics of 
the annihilating quarks are quite different for the two cases. The final state dileptons are 
expected to carry information on the environment of the annihilating quark-antiquark pair 
and therefore serve as a probe of the thermodynamical state of the quark-gluon plasma. 

The rate for the production of dileptons with an invariant mass Q per unit four-volume 
in a thermalized quark-gluon plasma is sensitive to the temperature T of the system and 
can be written in the form p2|,p^ 



(4.1) 



dQ'd'^ ' %=6:-A ^ ^ 2(271)4^^ ^ g2 ^ i^T. 
where crq{Q) is the lowest order qq — >• cross section at the center-of-mass energy Q given 
by i 
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m;, rriq are the rest masses of the lepton / and the quark q respectively, and Ki is the 
modified Bessel function of first order. 

In order to be meaningful experimental signals for the detection of QGP states, the 
thermal dileptons must be clearly delineated from other sources of dileptons, particularly 
the Drell-Yan background. It becomes important therefore to have good estimates of the 
dilepton production from the different sources at matching levels of accuracy. It is hence 
necessary to modify equation to incorporate the higher-order QCD corrections into 

the rate for dilepton production in the quark-gluon plasma. 

The interpolation from low to high relative velocities extends the usefulness of the cor- 
rection factors Eq. (|2.20|) . Therefore, we may introduce the same correction factor K^^^ as 
used in the DY case to account for high order QCD effects. The modified rate of dilepton 
production becomes 

The virtual photon annihilation mode selects the annihilating qq pair to be in the color- 
singlet states and thus the interaction is attractive. The distortion effect leads to an en- 
hancement factor K^'^^ as shown in Fig. 1. Because the coupling constant increases and the 
relative velocity decreases as the invariant mass decreases, the correction factor rises consid- 
erably with the decrease of invariant mass. The effect of the higher order QCD corrections 
is to enhance the tree-level dilepton cross section by a factor of about 4 at ^/s = 1 GeV, 
and by a factor of about 3 at ^/s = 2 — 3 GeV. Higher order QCD corrections will increase 
the strength of the dilepton signal and may serve to enhance the prospects for its detection 



in ongoing and planned experiments [33 



It may be remarked that in the literature, a correction K factor has sometimes been used 
for thermal production as [1 + (as/7r)(l + aT^/Q"^)] (Eq. (5.1) of Ref. [P^]), which coincides 



approximately with K^-l"^ that includes the first-order correction for the reverse qq 
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process pj, but is smaller than the correction factor K^'^\ Recent calculations, for the first 
order correction, using the thermal mass, predict a much larger K factor similar to our 
in magnitude p2 |. 

The question of screening in the quark gluon plasma requires a comment at this stage. In 
the plasma, the color charge of the constituents of the plasma is subject to Debye screening 
which is characterized by the Debye screening length A/? depending on the temperature 
T. For the quark-gluon plasma with a flavor number Nf and A^^^ = 3, the lowest-order 



perturbative QCD theory gives 



A^(PQCD) = ^ _ . (4.4) 

3 



For a coupling constant = (7 /47r = 0.5 and Nf = 3, the Debye screening length at a 
temperature of 200 MeV is about A^ ~ 0.4 fm. On the other hand, the electromagnetic 
annihilation into dileptons occurs within an extremely collapsed space zone characterized 
by the linear q-q distance ~ a/^fs which should be compared with the the Debye screening 
radius Xd in the plasma. As the annihilation distance is about 0.0029 fm for the 

annihilation of a light quark- ant iquark pair at 0.5 GeV, and is much smaller than the Debye 
screening length A/?, the interaction between the quark and the antiquark is not much 
affected by Debye screening in the region where annihilation occurs. Therefore, we expect 
that our correction factor for dilepton production by qq annihilation in the plasma will not 
be modified much by the addition of Debye screening corrections. 



V. HEAVY-QUARK PRODUCTION 

Heavy quark production in nuclear collisions, whether routed through intermediate QGP 
states or not, occurs either by qq annihilation (gg QQ) or by gluon-gluon fusion, {^gg — »• 
QQ). In the first case, the final state heavy quark-antiquark pair are required to be in a 
color-octet outgoing state, due to the color-octet nature of the intermediate gluon states, 
while in the latter, they may emerge in either the color-singlet or color-octet states. We 
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shall limit our investigation in this paper to the former process, and discuss the gluon-gluon 
process in a future publication. 

The cross section for qq QQ can be written at tree-level as [jl| 

In the qq —>■ g* QQ process, there are two strong-interaction vertices, qqg* and g*QQ. 
Both of these require QCD corrections, along with the other radiative corrections involving 
the external quark lines and their connections with the virtual gluon propagator. 

In the present work, we attempt to correct for the two interaction vertices simultaneously, 
to include the wave-function distortion effect at low relative velocities and additional real 
gluon radiation at high relative velocities, in the same way as in our earlier corrections for 
the DY process. 

The effective interaction between the quark and the antiquark can be written for each 
vertex as in section II, and we can similarly write a K factor for each vertex. The first 
vertex, qqg*, occurs with the annihilation of a qq pair in the initial state and is associated 
with the K^^\q) correction factor. The vertex is identical to the DY case except for the 
replacement of the virtual photon by the virtual gluon. Because of the color-octet nature of 
the intermediate gluon, the quark-antiquark pair must likewise be in the color-octet state. 
The correction factor K^^'^ (g) for the annihilation of this color-octet state corresponds to an 
effective repulsive potential with the color factor Cf equal to —1/6. As we note from Fig. 
1, it has the value of about 0.9 at high energies. 

The second vertex, g*QQ, occurs with the production of the QQ pair in the final state and 
is associated with the K^^^Q) correction factor given by Eq. ( p.20| ), having the appropriate 
function f^^\v) of Eq. (|2.20|) . Since the color factor is same as the first vertex, this correction 
is also suppressive. 

Corrected for wave function distortion at both vertices, the cross section ( |5.1| ) is amended 

to 

a{s)=K^''\q)K^f\Q)^[l-^j ^1-— Q(^l + 2 ' ^ Q + j . (5.2) 
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In the range of invariant mass close to the masses of the heavy quarks in question, the 
ideal analysis envisages correct inclusion of the heavy quark masses. The vertex factor 
K^^\Q) accordingly acquires a sensitivity to the mass. In Fig. 2, we display the correction 
factor K'^f'> for the production of qq pairs of various flavors in color-singlet or color-octet 
states, calculated with A = 0.3 GeV. The values of the correction factors are constants at 
relativistic relative velocities, but begin to depart from the constant values as the invariant 
mass approaches the threshold, corresponding to lower relative velocities between the quark 
and the antiquark. Near the threshold, the enhancement of the color-singlet states due to 
distortion is much more substantial than the corresponding suppression of the production 
of color-octet states. 

How does the Debye screening affect the correction factor? The annihilation or pro- 
duction cross section in the qq — > QQ process involves the strong coupling constant and is 
associated with an annihilation or production length of the order of as/i/i (c/. Eq. (|5.1| )) 
which should be compared with the Debye screening length A/? of about 0.4 fm at T = 200 
MeV. The effect of Debye screening is important when as/\/s » Ad; it is unimportant 
when as/\/s « Xo- Accordingly, for quark-antiquark production involving the annihila- 
tion or the production of light qq pairs through a virtual gluon intermediate state, the effect 
of Debye screening is important when ^/s « 0.25 GeV, assuming a coupling constant of 
Us = 0.5. It is unimportant when ^/s >> 0.25 GeV. Thus, the correction factors can be 
used for qq pairs with energies y/s much greater than 0.25 GeV. 

VI. CONCLUSIONS AND DISCUSSIONS 

At low energies of relative motion, the quark-antiquark relative wave function becomes 
distorted in the field of the color potential acting between a quark and its antiquark part- 
ner The distortion is strong in the vicinity of their zero separation and leads to a 
significant modification of the cross section for their annihilation or production. When the 
quark-antiquark pair is in its color-singlet state, the interaction is highly attractive, and 
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the cross section is much enhanced. On the other hand, when the quark-antiquark pair 
is in its color-octet state, the interaction is repulsive and the cross section is suppressed. 
This enhancement or suppression can be quantified in terms of a vertex correction factor K. 
The distortion effect does not depend on whether the quark-antiquark pair is annihilated or 
produced. At high energies of relative motion such that the quark masses can be neglected, 
the addition of real soft gluon emission leads to well-known PQCD vertex correction fac- 
tors which are different for qq annihilation or production Following the interpolation 
procedure of Schwinger to connect the vertex correction factor for low relative velocities, 
(obtained from wave function distortion), to the vertex correction factor for high relative 
velocities, (obtained from perturbative QCD), we can interpolate the correction factor to 
make it applicable to the whole range of relative velocities. The correction factor enables 
us to correct tree-level QCD calculations, when a qq pair is annihilated or produced. Its 
inclusion improves both our understanding of the underlying physics and the accuracy of 
the reaction cross sections. 

We wish to emphasize here the usefulness of our unified description to correct the lowest- 
order results in different environments. As our K^^'^^ does not dependent on the structure 
functions or constituent distributions in qq annihilation or the mode of qq production, it does 
not have to be recalculated for every individual environment. It is sufficient to ensure that 
the process considered is described by basic diagrams involving qq annihilation or production. 

As a test of the reliability of the correction factor Eq. ( p.20|) , we apply the correction 
factor for qq production to study the ratio R = (T(e+e^ hadrons)/cr(e+e^ n^ji^). 
For heavy quark production in the e+e~ — > 7* — > QQ process, Q and Q are produced in 
the color-singlet final state. The tree-level cross section must be multiplied by K^^^ of Fig. 
2 to take into account higher-order QCD corrections. The results of Fig. 2 shows a very 
large enhancement near the threshold of heavy-quark production. This large enhancement 
is indeed observed and the results Fig. 2 give good agreement with experimental ratio R = 
a{e^e^ hadrons)/cr(e+e^ fi^fi^) [§§]. The good agreement lends support to the 
usefulness of Eq. ( p.20|) for its application to other reaction processes. 
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Accordingly, one can use the correction factors to improve the tree-level estimates of 
many cross sections. For example, in dilepton production via the Drell-Yan process in 
nucleon-nucleon collisions, the correction factor K^^^ is about 2.2 from Fig. 1 at high en- 
ergies which agrees with experimental values of 1.6 to 2.8 In the quark-gluon plasma, 
dilepton production from the collision of quarks and antiquarks in the plasma occurs through 
the reactions uu and dd l^l^ . Because of our interpolation using Schwinger's 

procedure, the semi-empirical correction factors of Eq. ( p. 201 ) is applicable to these pro- 



cesses. The tree-level cross sections for uu l^l^ and dd — ^ l^l^ should be corrected by 
multiplying them with the correction factor i^^*^ , to take into account the initial-state u-u or 
d-d color interactions and radiations. The results of Fig. 1 indicates that this multiplicative 
correction factor for dilepton production in the quark-gluon plasma is significantly large. 
It increases as the dilepton energy decreases. The action of of the correction factor is to 
enhance the tree-level dilepton cross section by a factor of about 4 at ^/s = 1 GeV, and by 
a factor of about 3 at ^/s = 2 — 3 GeV. 

In the process of qq —>■ ss in the quark-gluon plasma, where g is a light u or d quark, 
the tree-level cross section should be multiplied by the correction factors K^'^\qq)K^^\ss) 
to take into account the initial- and final-state interactions and gluon radiations. Because 
these reactions take place when the quark and antiquark are in the color-octet state, the 
interaction is repulsive, and the two correction factors are less than unity. The deviation 
from unity is great near the ss threshold. Another application of the present investigation is 
the use of the correction factor to provide a better estimate of the production of top quarks 
by the e~^e~ annihilation, which will be the subject of a separate publication [^. 

In the plasma, the color charge of the constituents of the plasma is subject to Debye 
screening. We shall investigate quantitatively how the correction factors may be modified by 
the Debye screening. One can get a qualitative understanding whether the Debye screening 
may lead to a large modification on the correction factor of Eq. ( p.20| ) by comparing the 
magnitudes of length scale, over which annihilation or production takes place, with the 
Debye screening length. If the length scale for annihilation or production is much greater 
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than the Debye screening length, then the color charge of one interacting particle is effectively 
shielded from the other particle, when the quark and the antiquark come to the region of 
annihilation or production. In this case, the conrrection factor of Eq. ( p. 201 ) calculated with 



no Debye screening will be much modified when the Debye screening is taken into account. 
On the other hand, when the length scale for annihilation is much smaller than the Debye 
screening length, as in the case of dilepton production in the quark-gluon plasma, the effect 
of Debye screening is small. Accordingly, the Debye screening will not significantly modify 
the correction factor of Eq. (|2.20|) for dilepton production, and for the reaction qq QQ 
with ^/s >> 0.25 GeV. The Debye screening is important for the reaction qq — > QQ with 
^/s << 0.25 GeV. Therefore, it remains reasonable to apply the present results of the 
correction factor without considering Debye screening to a substantial domain of reactions 
in the quark-gluon plasma. 

An important part of the dynamics in the plasma or partons involves the reaction with 
gluons. The evaluation of the complete set of reaction cross sections in the plasma or partons 
will require the investigation of next-to-lowest order corrections to gluon reaction processes, 
which will be the subject of our future studies. Just as for a qq pair, initial- and final-state 
gluon-gluon and gluon-quark interactions are expected to lead to large corrections to cross 
sections for reactions involving gg or gq pairs. A gg pair can form many different color 
multiplets and the interaction between a gluon and another gluon is attractive when they 
are in their color-singlet and color-octet states. It has been argued that gluon dynamics may 
be described as massive spin-1 fields with the mass generated dynamically through strong 



gluon- binding forces [36|. Investigations on the gluon-gluon and gluon-quark correction 



factors using a gluon exchange potentials as in [p6| , p7| will be of great interest. 

Future refinement of the present investigation can be carried out by including the vari- 
ation of the coupling constant in modifying the spatial dependence of the color potential 
between a quark and its antiquark partner in the plasma. The correction factor we have 
obtained is based on an approximate simplified description by treating the two-body prob- 
lem as a single-particle problem with the proper relativistic kinematics. A better two-body 
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Dirac equation based on Crater and van Alstine's constraint dynamics p6[ has already been 
worked out and used to examine electromagnetic interactions Its application to the 

present problem involving a quark and an antiquark in a Coloumb-like interaction will be 
of great interest. Furthermore, the solution of Eq. (6) is actually obtained by neglecting the 
0(agg/r^) term in the second-order Schrodinger eqaution derived from the single-particle 



Dirac equation [241. It will be useful to include this term in future investigations. 
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FIGURES 

FIG. 1. The vertex correction factor for the annihilation of a qq pair of various flavors. 
The upper three curves are for color-singlet states and the lower three curves for color-octet states. 

FIG. 2. The vertex correction factor K^^^ for the production of a qq pair of various flavors. 
The upper three curves are for color-singlet state and the lower three curves for color-octet states. 



25 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/94 1 2349v 1 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/94 1 2349v 1 



Correction factor for qq production 




Fig. 2 



Correction factor for qq annihilation 



10.0 










9.0 

ft n 




t Ml M 

UU 

1 dd 


1 

1 cc 

1 


- 
- 
- 


7 n 
/ .u 




\ss 




- 


6.0 






1 




5.0 










4.0 






\ 
\ 

\ 

^^^^^^:=:.:i.,:,__ (^olor-singlet 


- 


3.0 






- 


2.0 










1.0 


uu 


dd ss 


color-octet 




0.0 


r 




cc 

1,1,1 





0.0 2.0 4.0 6.0 8.0 10.0 



(sr (GeV) 



